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ABSTRACT 

We present the AGILE gamma-ray observations in the energy range 50 MeV - 10 GeV of the super- 
nova remnant (SNR) W44, one of the most interesting systems for studying cosmic-ray production. 
W44 is an intermediate- age SNR (~ 20, 000 years) and its ejecta expand in a dense medium as shown 
by a prominent radio shell, nearby molecular clouds, and bright [SII] emitting regions. We extend our 
gamma-ray analysis to energies substantially lower than previous measurements which could not con- 
clusively establish the nature of the radiation. We find that gamma-ray emission matches remarkably 
well both the position and shape of the inner SNR shocked plasma. Furthermore, the gamma-ray 
spectrum shows a prominent peak near 1 GeV with a clear decrement at energies below a few hun- 
dreds of MeV as expected from neutral pion decay. Here we demonstrate that: (1) hadron-dominated 
models are consistent with all W44 multiwavelength constraints derived from radio, optical, X-ray, 
and gamma-ray observations; (2) ad hoc lepton-dominated models fail to explain simultaneously the 
well-constrained gamma-ray and radio spectra, and require a circumstellar density much larger than 
the value derived from observations; (3) the hadron energy spectrum is well described by a power-law 
(with index s = 3.0 ± 0.1) and a low-energy cut-off at E c = 6 ± 1 GeV. Direct evidence for pion 
emission is then established in an SNR for the first time. 

Subject headings: ISM: supernova remnants (W44) — cosmic rays — acceleration of particles — 
gamma rays: general 



1. INTRODUCTION 

Providing an unambiguous proof of the cosmic-ray ori- 
gin until now has been elusive, despite m any decades 



of attempts and controversial claims ( e.g., Fermi 1949 
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osmic-rays are mainly 
protons and heavy ions (hadrons) and, in a few percent, 
electrons and positrons. Supernova remnants (SNRs) are 
ideal candidates for the cosmic-ray production up to en- 
ergies near E^nee = 10 15 eV. The SNR energy output in 
the Galaxy can indeed supply the energy budget neces- 
sary to maintain the present population of cosmic-rays. 
Furthermore the observations of ultra-relativistic elec- 
trons support the hypothesis that also protons are accel- 
erate d in these objects (for a recent review see |Reynolds| 
2008, and references therein). Proving the fact that the 
SINK origin of hadronic cosmic-rays is difficult because 
of the complexity of the SNR-environment interaction. 
From an observational point of view, a direct proof can 
be given by an unambiguous detection of the gamma-ray 
emission expected from neutral pion decay in hadronic 
interactions. However radiation from co-spat ially ac- 
celerated electrons can mask and sometimes overcome 
the expected neutral pion decay signature of proton/ion 
emission in the 100 MeV-a few TeV energy range. Re- 
cent analysis suggests that several gamma-ray observa- 
tions of SNR s can be understood in terms of acceler- 
ated hadrons dCastro & Slane 2010; Giuliani et al.||2010[ 
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Abdo et al.|201Ub|c|al|2UU9[ ). However it is currently not 



possible to exclude that the observed 7-ray emission is 
produced by leptons alone. The discrepancies between 
leptonic and hadronic models are expected to be more 
evident at low energies (50 - 100 MeV). Gamma-ray as- 
tronomy in this energy band is very challenging because 
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of high background-noise flux and of the strong multi- 
ple scattering suffered by electrons originating from 7- 
ray events. The AGILE/GRID instrument (calibrated 
in the 50 MeV - 10 GeV band), however, has already 
shown its ability to provide a n energy spectrum starting 
at 50 MeV for bright objects ([V ercellon eet al.|2009| |Giu-| 
liani et aT]|2010 Vittorini et al.||201l|) . In this paper we 
report on a low-energy 7-ray and multiwavelength spec- 
trum for the SNR W44 in order to constraint the emitting 
particle spectrum and discriminate between leptonic and 
hadronic models. 

2. THE SUPERNOVA REMNANT W44 

SNR W44 (G34.7 -0.4) is a well studied middle-aged 
(~20,000 yr) SNR located in t he Galactic disk at a dis- 
tance of - 3 kpc from Earth flClark et aL]|1976| |Wol- 



et al. 



1991[ ) . W44 is an ideal system to test the 
presence of accelerated hadrons and the in terplay be- 
tween hadro nic and leptonic models. Radio (Castelletti 
et al.||2007[ and references therein) and X-ray (|\Vat-| 



son et al. p.983 ) mapping of the SNR show a roughly 
elliptical shocked shell and a centrally peaked emission 
respectively. In the IR band ( |Reach et al.|2005| ) it is evi- 
dent that the shell is expan ding into a dense surroun ding 
medium (n ~ 100 cm -3 ). |Wolszczan et al. (1991) dis- 
covered the radio pulsar PSR B1853 +01 with distance 



and age c ompatible w ith the SNR. [Wootten et al.| ( |1977 ) 
and then |Rho et aT] ( |1994| ) found a massive molecular 
cloud interacting with the south-eastern side of the rem- 
nant. Evidence for a more complex system of massive 
MCs and for their interactions with the remnant, shown 
by some features typical of a s trong shock, was giv en by 
|Seta et al.J ( |2004[ ) and then by |Reach et~at| ( |2005[ ). The 
MC-SJNR interactions were confir med by the maser OH 
(1720 MHz) emission reported by Claussen et al. (1997) 
and then by jHoffman et al. ( 2005 ) . 

The first estimation of the spectral radio index varia- 
tio ns as a function of positi on over the remnant was done 
by |Castelletti et aT] ( |2007| ): the eastern limb spectrum 
was consistent with a diffusive shock acceleration model 
and the spectrum flattening in the westernmost arc con- 
firmed the MC-SNR interaction observed in IR and op- 
tical band. 

Gamma-ray emission from this SNR has been detected 
by the Fermi /L AT instrument at energy E>200 MeV 
flAbdo et al.||2010a| , suggesting the presence of acceler- 
ated protons interacting with the surrounding medium. 

3. DATA ANALYSIS 

AGILE-GRID data were analyzed using the AGILE 
Standard Analysis Pipeline. We used 7-ray events fil- 
tered by means of the FM3.119-2 AGIL E Filter Pipeline 
(as described in |Vercellone et aL] |2008). In order to dis- 
criminate between background events and gamma rays, 
the GRID and anticoincidence system (ACS) signals are 
processed , reconstructed and selected by a dedicated 
software ([Giuliani et al. 



[20061). W e used the most re- 
cent versions of the diffusion model (|Giuliani et al.| 2004) 
and of the calibration files, available at the ASDC site 
(www.asdc.asi.it). We created counts, exposure and 
Galactic background gamma-ray maps with a bin-size 
of 0.02° x 0.02°. In order to derive the source average 
flux and spectrum we ran the AGILE point so urce anal- 
ysis software ALIKE ( |Bulgarelli et al. 2011 based on 



the maximu m likelihood technique described in Mattox 
et al. (1993)) over the whole observing period 2007 July 
- 2011 April. Both statistic and systematic uncertainties 
are taken into account. The spectrum was obtained by 
computing the 7-ray flux in six energy bins selected with 
the aim to have a significance a > 4. In order to study 
the source morphology, we obtained an intensity map in- 
tegrated over the energy range where the GRID angular 
resolution is optimal (E>400 MeV). 

4. RESULTS AND DISCUSSION 

AGILE detects SNR W44 with a significance of 15.8 a 
as an extended source. Figure [T^i shows AGILE gamma- 
ray intensity map above 400 MeV of the W44 region 
with the 324 MHz VLA radio contours. The gamma- 
ray morphology remarkably resembles the quasi-elliptical 
pattern of the interior of the radio shell, especially coin- 
ciding with the radio brightness enhancements toward 
the north-west and south-east r egions. Both the radio 
pulsar PSR B1803+01 position (|Petre et al. 2002) and 
its (small) pulsar wind nebula (Giacani et al 1997 5 and 
references therein) are inconsistent with the gamma-ray 



morpho logy detected by AGILE. Moreover, |Abdo et aT 
( |2010a[ ) excluded the presence of a pulsation in the 7-ray 
signal. Figure shows the CO emission at a kinematic 
velocity compatible with the distance of W44, tracing 
the presence of MCs in the W44 surroundings, together 
with the gamma-ray contour levels. It can be inferred 
that the south-eastern side of the 7-ray source overlaps 
with the MC - SNR interaction region. In the northern 
part of the shell, the 7-ray and CO emissions are not 
correlated, however many studies of the surrounding in- 
terstellar medi um showed the pr esence of dense gas not 
traced by CO ( |Reach et al.||2005 ). Figure [1J displays an 
SII map overlaid with X-ray and gamma-ray contours. 
The strong sulfur [SII] emission, along with Ha emission, 
indicates the presence of shocked gas ( |Draine fc McKee] 
1993) ' 
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4.1. The gamma-ray spectrum 

Figure [2] shows the AGILE W44 photon energy spec- 
trum for the whole range 50 MeV- 10 GeV. The measured 
flux above 400 MeV is F = (16.0 ± 1.2) x 10~ 8 pho- 
tons cm -2 s -1 . In this figure are a lso shown Fermi/LAT 
spectral points ( Abdo et al.pO lOa) over an energy range 
0.2-30 GeV. In the band where the spectra overlap, the 
sets of data are compatible at la. In this regard it is im- 
portant to stress how AGILE is able to detect the emis- 
sion from W44 in the energy range 50 MeV-300 MeV 
extending the spectrum to energies substantially lower 
than those previously obtained. This spectrum shows a 
clear decrement at photon energies lower than 400 MeV, 
confirming the expectations based on neutral pion emis- 
sion from accelerated protons/ions: a peak energy near 
1 GeV for hadron energy spectra flatter than E~ 2 (e.g., 
Aharonian 2004). The gamma-ray spectrum becomes 
steep at higher energies with a photon power-law index 
a ~ 3±0.1, in ag reement with previous measurements 
( |Abdo et al.|2010a[ ). 



18 As pointed out by |Rho et al.||l994| in their X-ray and optical 
study of W44, the optical filaments and the X-ray image showing 
locally bright emission clumps along the filaments suggest that 
both are produced by the interaction between the SN shock front 
and regions of enhanced ambient density. 
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Fig. 1. — (a): AGILE gamma-ray intensity map (in Galactic coordinates) of the W44 region (1.1 x 0.75 degrees) in the energy range 
400 MeV - 3 GeV obtained by integrating all available data collected during the period 2007 July and 2011 April. The color bar scale is 
in units of 10 — 5 photons cm -2 s — 1 pixel -1 . Pixel size is 0.02 degrees with a six-bin Gaussian smoothing. Green contours show the 324 
MHz radio continuum flux density detected by the Very Large Array. The white cross indicates the position of the pulsar PSR B1803+01. 
Source A does not appear to be associated with W44. (b): combined CO data from the NANTEN Observatory superimposed with the 
AGILE gamma-ray data contours above 400 MeV (in white) of the W44 region (map in Galactic coordinates). CO data have been selected 
in the velocity range 35 - 45 km s — 1 corresponding to a kinematic distance compatible with the W44 distance, (c): the [SII] image of SNR 
W44 obtained in 1993 with the 0.6 m Burrell-Schmidt telesc ope at Kitt Pe ak National Observatory with a matched red continuum image 
subtracted to emphasize the faint nebulosity (adapted from [Giacani et al||1997| >. The overlaid contours trace the gamma-rays emission 
detected by AGILE between 400 MeV and 3 GeV (in red) and the A-ray emission detected by ROSAT in the 0.1 - 2.4 keV range (in blue). 



4.2. Lepton- dominated models of emission 

Leptonic-only models of gamma-ray production have 
to satisfy the very well determined spatial and spectral 
constraints provided by the radio, optical, and gamma- 
ray emissions. For both the Bremsstrahlung and in- 
verse Compton cases we tested the class of leptonic-only 
models attempting to reproduce the observed gamma- 
ray spectrum for different input leptonic spectra (see 
Table 1). We used a phenomenological approach 
based on two evidences: a) the synchrotron spectrum 
implies that the radio electron distribution is well de- 
scribed by a power-law over a wide range of energies, 
b) the discontinuity in the slope of the gamma-rays 
spectrum implies a discontinuity in the emitting par- 
ticle spectrum. We assumed three electron distri- 
butions that can, in principle, describe this behavior: 
1) a power-law with a high energy cut-off, F\e{E) = 



K e E- 



Hendrick et al.|2001 ), 2) a power-law 
with a low energy c ut-off, F 2 e(E) = K e E~ p e~^ , (as in 



: , (as in 



Gabici et al.||2QQ9[ ), 3) a broken pow er-law, F 3 e(E) = 



Zirakashvili et al. 



2007), where E c is the cut-off energy and K e is the nor- 
malization constant. 

Since the distributions with a cut-off (1 and 2) fail to 
reproduce simultaneously both the radio and 7-ray spec- 
trum, we refer to the distribution (3). The best fit to 
the gamma-ray data is obtained with E c — 1 GeV, and 
indices p\ = and p 2 = 3.3 above and below E C1 re- 
spectively. Synchrotron emission originating from this 
distribution can be evaluated for different values of the 
average magnetic field. Figure [3] shows the case of the 
most "favorable" leptonic-only model characterized by 
B = 20 jiG (other cases turn out to be even less fa- 
vorable). We find that, at high frequencies, the cal- 
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Fig. 2. — Combined AGILE (red) and Fermi/LAT (green) spec- 
tral energy distribution (SED) for SNR W44. AGILE points are 
in the range 50 MeV- 10 GeV divided into six energy intervals. 
Fermi/LAT data span the energy range 0.2 - 30 GeV (from Abdo 
et al., 2010). 

culated synchrotron spectrum is in strong disagreement 
(factor larger than 4) with the ra dio emission produced 
co-spatially to the gamma-ray one ( [Castelletti et al.| 2007, 
and Figure]^). Furthermore, the inferred average density, 
n= 300 cm -3 , is too large (by a factor of three ) com- 
pared with the circumstellar medium constraints (Reach 
et al. 2005); a lower value for n would be incompatible 
with the gamma-ray spectrum fitting. Similar or even 
stronger contradictions with the multiwavelength data 
apply to other leptonic-only models that we systemati- 
cally explored for a large variety of parameters. 
In the case of inverse Compton dominated models, two 
sources of soft photons are available: the cosmic back- 
ground radiation (CBR) and the interstellar radiation 
field (ISRF). In the first case, a second peak in the 
gamma-ray spectrum is unavoidably expected with a 
peak energy E max ~ 1 TeV, in contradiction with the 
upper-limits obtained from TeV Cherenkov telescopes. 
In the case of interaction with ISRF, instead, the calcu- 
lated synchrotron peak is not compatible with the radio 
continuum data for any reasonable value of the magnetic 
field in the SNR shell. We can then reliably exclude 
leptonic-only models of emission for SNR W44. 

4.3. Hadron- dominated models of emission 

It is interesting to determine the main physical pa- 
rameters of an emission model dominated by hadrons in 
the gamma-ray energy range. The gamma-ray emission 
was derived assuming that protons interact with the nu- 
clei of the ambient medium thro ugh pp interactions and 
then radiate through tt° decay ( jKelner et al.||2006[ ). In 
order to fit the gamma-ray spectrum we tested, as pro- 
ton energy distribution, a power law, a broken power 
law and a power law with cut-offs. The distribution pro- 
viding the best-fit of the gamma-ray spectrum turned 
out to be a power-law with a relatively steep spectral 
index and a low-energy cut-off, F P (E) ~ E~ p e~ Ec / E , 
with p = 3.0 ± 0.1 and E c = 6 ± 1 GeV. The inferred 
total energy of hadrons producing the observed gamma- 
ray spectrum is Eh = 10 49 — 10 50 erg (depending on the 
average magnetic field and local densities, see Table 1) 



corresponding to a fraction of the total SNR shell kinetic 
energy e ~ 0.01 — 0.1. This value is a lower limit on 
the total energy of hadrons accelerated by W44 during 
its entire life, because the most energetic hadrons, likely 
accelerated during the early epochs of the SNR life, can 
escape from the accelerat ion site on timescales shorter 
than the ap e of this SNR (Ber ezinskii et al^]|1990| Gabici 
et al.||2009 ). Electrons produce the radio continuum and 
shell-like features by synchrotron emission in an aver- 
age magnetic field in the range B ~ 10 — 100 fiG. They 
also produce Bremsstrahlung and inverse Compton com- 
ponents that are subdominant in the gamma-ray energy 
range for a standard value of the electron/proton number 
ratio x ~ 0.01. Our best "hadronic model" is shown in 
Figure [4j We note that the spectral index found in this 
analysis is quite steep in comparison with the expecta- 
tions for_Jth^_spe£trum of particles accelerated in SNRs 
(e.g. Reynolds 2008) and that the cosmic-ray spectrum 
can be altered (softened) by the interaction with sites 
where the target g as is confined, for example by molecu - 
lar clouds (see e.g. Gabici et al.|2010 Ohira et aLj|2011 ). 
A comprehensive theoretical interpretation of this fact is 
beyond the scope of this work and will be addressed in a 
forthcoming paper. 
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Fig. 3. — Most favorable of the leptonic-only models, character- 
ized by B=20 /iG and n= 300 cm -3 . The green curves show the 
electron synchrotron (dashed curve), Bremsstrahlung (solid curve), 
and IC (dotted curve) contributions. This model is in contradic- 
tion with the radio data, and requires an average value of the gas 
number density that is too large compared with the value deduced 
by radio and optical observations (n ~ 100 cm -3 ). 



5. CONCLUSIONS 

W44 turns out to be an extremely interesting SNR 
whose environment and shocked material configuration 
favor a detailed testing of hadronic versus leptonic emis- 
sion. We find that models of accelerated protons/ions 
interacting with nearby dense gas successfully explain 
all the observed multifrequency properties of W44. The 
total number of accelerated hadrons in W44 is consistent 
with an efficiency of a few percent in the conversion of 
SNR kinetic energy into cosmic-ray energy. 

Within the hadronic scenario, the broad-band spec- 
trum is successfully modeled by the gamma-ray emission 
resulting from the decay of neutral pions produced by ac- 
celerated protons/ions with a spectral index near 3 and a 
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Fig. 4. — Theoretical modeling of the broad-band spectrum of 
SNR W44 superimposed with the radio (data points in red color) 
and gamma-ray data of Figure[2](in blue color). TeV upper limits 
are also shown.. This is an hadronic model characterized by B 
=70 fiG and n =100 cm -3 . The yellow curve shows the neutral 
pion emission from the accelerated proton distribution discussed 
in the text. The green curves show the electron contribution by 
synchrotron (dashed curve), Bremsstrahlung (solid curve), and IC 
(dotted curve) emissions. The red curve shows the total gamma- 
ray emission. 



cut-off energy E c ~ 6 GeV (Figures [5] and [4]). Electrons 
(a few percent in number) contribute to the synchrotron 
radio emission and to a weak Bremsstrahlung compo- 
nent. It is interesting to note that the spectral cut-off 
energy E c and the steep power-law index may indicate 
either a "diffusion" effect of high-energy hadrons, or a 
suppression of efficient particle acceleration in dense en- 
vironments (Uchiyama et al. 2010 Malkov et al. 2011) 
as expected in W44. Our results on W44 are comple- 
mentary with those obtained in other SNR s, most no- 
ta bly IC 443 flTavani et aI][20Tol) and W28 flGiuliani et| 
al. 2010). Also in these SINKs the complex interaction 
oi cosmic-rays with their dense surroundings produces 
gamma-ray spectra with high-energy cut-offs observed 
in the range 1-10 GeV. 

We acknowledge stimulating discussions with F. Aha- 
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Space Agency (ASI) with scientific and programmatic 
participation by the Italian Institute of Astrophysics 
(INAF) and the Italian Institute of Nuclear Physics 
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the ASI grant no. 1/042/10/0. G.D. and G.C. are mem- 
bers of CIC-CONICET (Argentina) and are supported 
by CONICET, ANPCyT and UBACYT grants. 
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Table 1: Electronic component of hadronic and "ad hoc" leptonic models 



Model 


< n > 
(cm- 3 ) 


< B > 


E h 
(erg) 


E e 
(erg) 


Pi 


P2 


Comments 


Hadronic 


20 


20 


2.1 x 10 50 


2.0 x 10 49 


1.7 


3.0 


• good agreement with data 

• n small 

• E h large 


100 


70 


3.3 x 10 49 


2.8 x 10 48 


1.7 


3.0 


good agreement with all data 


Leptonic 


300 


20 




7.3 x 10 48 





3.3 


• incompatible with the spectrum 

• too large medium density 

• too large E e 


5000 


200 




4.4 x 10 47 





3.3 


• incompatible with the radio spectrum: wrong 
slope, low- frequency excess 

• too large medium density 

• too large average B 

• electron cooling (Bremsstrahlung) timescale too 
short with respect to the SNR age 



